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Abstract 

We look at the Equivalence Theorem as a statement about the 
absence of polynomial infrared divergences when niw — > 0. We prove 
their absence in a truncated toy model and conjecture that, if they 
exist at all, they are due to couplings between light particles. 
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I Introduction 



The Equivalence Theorem [jT|, from a practical point of view, is a tool to 
simplify calculations in spontanously broken gauge theories. It states that, 
in the limit oo, amplitudes with external longitudinal W^s and Z^s 

can be replaced by the corresponding (unphysical) Goldstone boson ampli- 
tudes. With this replacement one gains a simplification in Lorentz indeces 
and calculations are simplified by the lack of gauge cancellations. 

The early incomplete proofs of the Equivalence Theorem within the stan- 
dard model were improved using power counting arguments 0, see 
also 0. One distinguishes two regimes, the light-Higgs regime with ^/s » 
rriH, TTiiv, nif and a heavy-Higgs regime with rriH ~ ^/i >> rriw, m/; in both 
cases excluding situations with exceptional momenta {{pi — pjY « E^) 
which could upset the power counting. The power counting proof of the 
Equivalence Theorem as given in works to all orders; it also shows that 
in the heavy-Higgs regime all Feynman diagrams with lines in loops other 
than Goldstone scalar are subleading, allowing for a consistent truncation of 
Feynman amplitudes to pure Goldstone dynamics. 

In this paper we discuss a possible flaw in the power counting proof which 
can occur if there are power-like infrared divergences. We argue that these 
do not occur in the minimal Higgs sector of the Standard Model, and use 
a toy model to argue that they are probably not present in its multi-Higgs 
extensions either. 

C. Grosse-Knetter's argument P] involves counting powers of mn and 
E = ^/s simultaniously. One starts with the well-known statement that 
< phys I TFi . . . Fn \ phys > = for the i?g-gauge fixing terms Fj = 

~ ' d^^i 0, connecting matrix elements of the Goldstone bosons 
04 to vector bosons with (unphysical) polarization oc Then one is to 
prove that the difference between the polarization vectors of these V^s and 
the longitudinal Vl vector bosons, = — ^ = ^^^^ ■ (1 | O)'^ gives 
only subleading contributions in the appropriate limit. A generic Feynman 
amplitude is, explicitely displaying the E and mu factors, a sum of terms, 
each in the form 

M=c-E^-^^ ■m^^'' - If, (1) 

where the constant c may depend on the my's, V$ is the added number of 
and $^ Higgs self-interactions, Ej is the number of external fermions and 
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is the number of v'^^s; the remaining part of the Feynman amphtude is in 
the general form 

^4L^ P---P /2) 

{qi - ml)... [q] - mj) ' 

where we collectively denoted by p and q the occurring linear combinations of 
external and loop momenta; m denotes both the heavy and the light massses. 

Eqn. ^ is a general formula for all processes and should include all internal 
lines present in the theory. As explained below, an assumption on its behavior 
in general at large scales can be used to prove the Equivalence Theorem. For 
this purpose, the relevant graphs are those which include at least one external 
longitudinal vector boson. 

If we now suppose that Ip is determined by the scale E and uih, we have 

lF = m^- /o(— ) + mg-i ■ mw ■ ) + • • • (3) 

with dimensionless fj (possibly containing logarithms of niw and of the renor- 
malization scale but no powers of ^— ) easy combinatorics shows that the 
total power of E and vtin-, counted simultaniously, in M. is at most 

iV = (2L + 2) - {Vd + 21^0 + V7 + E,) (4) 

where L is the number of loops, and Vj are the numbers of derivative 
and fermion couplings respectively; Vq = V — V.^, — Vd — Vf. This formula 
shows that the leading graphs are, at each loop level, those with E^ = (i.e. 
proves the Equivalence Theorem) and also shows that no vector or fermion 
lines are involved in loops in leading graphs. The former statement can also 
be proven true for the light-Higgs regime using a similar argument. 

A possible flaw in this argument can occur if negative powers of the light 
mass my enter the expansion (|^). That expression can be viewed as an 
infrared statement: up to possible logarithmic factors involving the renor- 
malization scale, we can fix the unit of dimension at E (or at mh (x E), to 
get a theory with my — >■ 0, wherein my acts as an infrared regulator. Then, 
the statement in (|^) has been transformed into the statement that the lead- 
ing part of a graph, fo, does not pick up polynomially divergent factors when 
my — 0. In many calculations of IR divergent quantities that use a small 
regulator mass (for examples, see P, p| among many others) one finds that 
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the divergences are only logarithmic in the regulator mass and this is the 
basis for the expansion in (|^). Any infrared divergences of the polynomial 
type in — would introduce additional factors of — and/or — and such 

^ my my I my 

terms would invalidate the above proof which is based on their absence. To 
see if any such factor actually upsets the Equivalence Theorem one needs 
to analyse by how many powers of my the graph in question is suppressed 
in (0). For exemaple, VlVl — > VlVl at one loop has (^-j 
as leading behavior. Any additional factor of ^ 



E or mw ) ' ^'^^ ^ ^ 

in, say, the box diagram 

of Fig. D, which itself is one of the leading graphs, would introduce a factor 
that breaks the Equivalence Theorem. 



my 




Wl 



Wl 



Figure 1: One of the leading graphs that might pick up IR divergences at 
1-loop level. 



A related problem could occur in a theory with several Higgses when some 
of them are heavy and at least one is light. Powers of — for the light Higgs 
could upset the power counting. 

What the argument in does is, in this light, that it shifts the question 
of 'breaking' the Equivalence Theorem to the question of the presence of 
severe (i.e. polynomial) IR divergences. 

This point of view is in accordance with the view that the Equivalence 
Theorem expresses the fact that Goldstone d.o.f 's are turned into longitudinal 
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vector bosons by a spontanously broken gauge transformation. With fixed 
E and mv — >■ we approach the point of phase transition where the t'Hooft 
gauge condition d^V^ = my/ ■ turns from a way to express with Vl into 
the physicahty condition on Vl. The m — hmit of a vector boson is a 
notoriously tricky problem. Whether the transition is smooth enough not to 
break the relations between amplitudes or not, will show up in the presence 
or lack of 'bad' IR divergences. 

In this paper we argue that such divergences do not occur in the pure 
Higgs sector, at least, for non-exceptional momenta. We try to relate the 
IR divergence found in scalar box graphs (see Eqn. (|5|) below to the Equiv- 
alence Theorem and find no such divergence when the external particles are 
massive. In the copnclusion we point out that a similar treatment is neces- 
sary for massless external particles which requires an analysis more similar 
to the case of the IR divergencesd in QCD. Experience with graphs involving 
massless particles of different spins shows that the infrared divergences of 
a diagram are generally determined by the denominators in the integrals of 
Eqn. (0) and the presence of spin affects only the numerators in that ex- 
pression. The result is that the general structure of the IR divergences is 
independent of the spin (see e.g. 0) and the reason is the following. IR 
divergences come from specific points in the loop integration region where 4- 
momenta of some massless propagators vanish. One can then expand the 
numerators around those points in k. We see that in terms of higher order 
in fc, IR divergences cancel and the remaining divergent part involves only a 
constant, i.e. a structure that is similar to the spin-0 case. This argument 
tells us that the absence of IR divergences in the scalar sector essentially 
proves the absence of divergences in the full model. One might hope that for 
some theories one could get special behavior from the numerators for special 
reasons, such as gauge invariance, that might remove an IR divergence that 
would otherwise be present. However, the scalar theory then would only be 
worse, so that by exploring the scalar theory one should be able to find the 
strongest divergences. 

The situation here is analogous to QED in that interactions of massive 
particles through massless photons (the latter now correspond to the W and 
Z) introduce only logarithmic IR divergences. To complete the proof of the 
Equivalence Theorem to all orders one should also prove that in a situation 
more reminescent of QCD, where massless gluons (now corresponding to W 
and Z) introduce IR divergences that are much harder to handle, polynomial 
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IR divergences are also absent - a question we do not address here. 

The form of IR divergences of Feynman graphs has been calculated in 
QED many times (see, e.g. Zwanziger 0), and is usually found to be loga- 
rithmic in the various IR regulators to any finite order. It has been calcu- 
lated PI for gravitons with similar results. However, it is well-known that in 
the nonrelativistic limit scalar box graphs, such as that in Fig. ^. pick up 
O terms for forward scattering 0: 

iMr^ i : (5) 

m?M{p + |m) ^ ^ 

clearly showing the type of behavior expected to break the Equivalence The- 
orem. It is worth to note, however, that these calculations use either a small 
off-shellness parameter Si = pi — ml or an explicit cutoff | p° | < A as 
an IR regulator, so they are not directly relevant to our case where a small 
regulator mass m should be used. 



M 



M 



M 



m 



m 



M 



M 



M 



Figure 2: One- loop graphs with powerlike — behavior. 

It is an interesting coincidence that the structure in Fig. |^. is not present 
in the Standard Model. The only heavy mass there is M — > mu and there are 
no HHW or HHZ couplings (the absence of the latter is a consequence of 
separate C and P-conservation in the purely bosonic part of the SM, where 
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J^^{Z) = 1 ). All other box graphs are less IR divergent because less 
interal lines can be put on shell in the same time. 

The obvious way to look for a similar divergence is in the two-doublet 
Higgs model which posesses one heavy Higgs M uih and at least one 
light 0^+ Higgs m < niz- Even though such an IR divergence is 

not a genuine because of the necessary lightness of ruh, it affects the 
amplitudes in the same way. 

In Sect. |T|we use Weinberg's method to analyse the 'truncated MSSM 
model' to all orders, dropping all particles from there except the heavy and 
the light O'^"*' Higgses, using m rrih as an infrared regulator. This model 
conserves the number of heavy i^'s, which so act as 'charged' particles while 
the light /I's act as uncharged scalar photons. The lack of the coupling be- 
tween light particles allows to extract the IR divergences in complete analogy 
to the QED case and we find only logarithmic divergences. The S'-matrix 
elements pick up factors 

S,. = S^;^ ■ exp 1^ ■ + ^:F) ■ log A| (6) 

where 2ng = Gmz is an MSSM coupling constant, m -C A <^ is an energy 
when we separate soft particles from hard ones. has contributions from 
pairs of incoming and from pairs of outgoing i/'s; Q has contributions from 
all pairs. As we will see in Sect. ^ ^ is canceled on the cross section level 
by real soft /i's; the T does not contribute to the cross section but it shifts 
the Coulomb phase by an ~ log — term, thus explaining why there is no 
Coulomb phase contribution from a pair formed of one incoming and one 



outgoing H. In Sect. |T| we calculate explicitely the IR divergent contribu- 
tions in HH HH to one-loop level and find complete agreement with the 
above conclusions. It is of some interest how these IR divergent terms are 
separated using our small-mass regularization: one needs a careful and long 
procerdure which is illustrated by describing the details for one particular di- 
agram. The mechanism of the cancellation of IR divergences, all logarithmic 
in m, between loop integrals and soft 'photons' attached to heavy external 
legs closely follows the corresponding mechanism in ee — > ee in QED, so it 
is straightforward to generalize the result of Sect. |I| for any other process 
in this toy model, for example, for hH — > hH or hh hh, analogous to 
67 — i> 67 or 77 77. 
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The total absence of powerlike divergences shows that the ^ divergence 
of the nonrelativistic forward amplitude is due to exceptional momenta. This 
statement is in compliance with that the coefficient of our log — divergence 
itself diverges when two particles are collinear (see Eqn. ([T^), 0). 

In Sec. m we show that in a simplified model with two scalars only no 
polynomial IR divergences occur. As a consistency check we show that the 
IR structure found is identical with the one necessary in order that soft 
Bremsstrahlung exactly cancels the IR divergences in loops. We then illus- 
trate this general argument in Sect. |T| with a calculation of the relevant 
1-loop amplitudes. In the conclusion we point out that the presence of cou- 
pled light particles (i.e. gauge boson self-coupling) introduces additional 
difficulities which require a separate investigation. 



II The 'truncated MSSM' model 

We use the H and h part of the MSSM as a toy model to illustrate why 
IR divergences are at most powers of logarithms in the small masses. This 
model has only one coupling (see Fig. ^ with a dimensionless G ~ gweak] 
the mz factor in the couphng should not be considered as suppressing the 
IR divergence: any deviation in the expansion of Ip from Eqn. (^ upsets 
the power counting. In this model, with M itlh and E kept constant and 
m rrih sent to zero, we calculate the IR divergences, closely following the 
argument in Weinberg [Q. 

All the IR divergences in an amplitude stem from a set of soft h exchanges 
between external legs of Feynman graphs. Their factorization [^] is due to 
the fact that these corrections are attached to all graphs in the same way. 
We have, for each attached soft exchange (see Fig. |[) 

Sp^ = Sfl-[l + j d'q.A{q)] (7) 

and 

= (2^7g2-l' + z.- 

i{i ■ 2ng) i{i ■ 2ng) 
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Figure 3: The only remaining vertex in the truncated MSSM. 

and T]j = +1 (or — 1) when the H is outgoing (incoming). Attaching all 
possible soft exchanges leads to a factorized amplitude 

S^^ = ■ exp J d'q ■ A{q) . (9) 

The integration is over a [—A, +A] range in with IR cut ^ A ^ i?; 
the resulting A dependence is compensated by a A-dependence in the hard 
h part Sf}^. In addition, by the usual \ q^ \> 2 cutoff we get, after some 
elementary integrations, with m = 0: 

^2 



9 VjVi 1 + A 

where f3ji is the relative velocity of two H^s 



^A = ^-E7^-f logf^ (10) 
8 ^ [pj ■ pi) 13 ji 1 - 



(5. 



\ KVj ■ Pir 



What we need is to use m > instead of A > 0; to calculate it we 
need a careful analysis of how complex singularities move around on the 
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Figure 4: Attaching a soft exchange to external hnes. 



plane. Closing the contour around them and integrating over angular vari- 
ables leaves us with one- dimensional integrals. Separating their IR divergent 
(i.e. when m — >• 0) parts is a tedious calculation; we only quote the result in 
the form of Eqn. with 



and (12) 
j^i with v,=m 

In a physical process, the exp \i ■ ^ ■ ^ ■ log ^} factor contributes to the 
Coulomb phase as we said in the introduction; the Q part goes into the cross 
section and gets compensated by real soft /I's. 

One may always add any number of indetectably soft /i's to the initial 
and/or final states as long as their total energy is less then the energy reso- 
lution of the measuring device. Although this has not much practical sense 
for nih ~ lO's of GeV^s, we still see that a Bloch-Nordsieck-type cancellation 
occurs. Attaching N soft /I's to the external legs, the amplitude picks up an 
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IR divergent factor 



N 



5'/3a(qi, • • • , (In) — S, 



Pa ■ 



HE 



i{i ■ 2'Kg) 



(13) 



These states should be added on the probabihty level as they represent 
orthogonal states. Using the factorization property as in ^ again, we get for 
the transition rate 



r/3a(qi, • • • , qAf) 



■E 



1 



N 2 

— ■ IT — — 

71=1 Qn 



(14) 



VjVi 



ll [{Pj ■ <ln) + Vj^^^] [iPl ■ Qn) + Vl 



Integrating this over the /I's' phase space with the restriction 



^ A 1 f+'^ da 

\n=l J J-co a 



sm cT ■ e A-z^n=i-^"- 
we have, with the total energy carried away by the /i' less than A: 



(15) 



Pa 



;< A) 

■ exp { 



S, 



Pa 



9_ 
16 Jo 





/■+°° da 


TT 


l-oo a 


dq ■ 









sm a- 



(16) 



E 



VjVi ■ ^^I^ 



4n 



(Pj ■ q) +Vj- 



{pi ■ q) + vr 



} 



The separation of IR divergences requires a hard mathematical procedure 
that we do not present here; in essence, we write the integral 

dq ■ q^ 



(17) 



dq ■ q^ 
E. 



(e^i^, _1) + 



dq ■ q^ 
E„ 



and prove that the first two terms do not contribute to the IR divergence. 
Our result is 



Pa 



Sf3a P - exp 



4(2^)^ 



G ■ log 



A 



m 



(18) 
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with the same Q as in Eqn. (JT^), proving the cancellation of all IR divergences 
in the transition probabilities. 

Expanding any of our results in the coupling constant G shows that, to all 
orders, the worst divergence is only a power of log ^. We note that although 
summing up all orders we certainly get a power like behavior, 

V = (19) 

we do not think this points to the breaking of the Equivalence Theorem; we 
are, after all, summing for a very particular set of diagrams. 



Ill A one-loop example 

In order to illustrate the results of Sect. ^ and also to see how the IR diver- 
gences of particular diagrams add up, we work out the infrared divergences of 
HH —>■ HH at one loop level and find complete agreement with Eqns. (|6|JT^). 
Fortunately, all graphs with IR divergences are UV finite, so we may ignore 
renormalization. It turns out that each individual divergent graph has a 
log — divergence and no cancellations occur. 

The IR divergent graphs are shown on Figs. In addition to these, 
one must include tree graphs and add all h mass insertions to them; these 
graphs turn out to be IR finite though. The calculation of each individual IR 
divergent part is too complicated to explain here; we briefly describe one of 
them (the one corresponding to Fig. |^ in the Appendix). We simply quote 
the result in terms of I graph- 

iMgraph = -71"^ ■ g'^-log—- Igraph (20) 



and 



m 



J(4b) = ~-^{u) + {t^u) (21.2) 
/(4c) = -- ■ <Pit) + (t ^ u) (21.3) 
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/(5a) 
I 



(56) = -- ■ <Pif) + {t^u) 



(21.4) 
(21.5) 



with E = Ia/s; (3 is the relative velocity /? 



1 _ 4Mi and 



log 



a/4M2 - X + 



^x{x - 4M2) Vil? 



(22) 




(a) 



/ \ 

^ X. 

(b) 




(c 



Figure 5: Infrared divergent box graphs in the truncated MSSM. 
The sum of all these terms gives (with spacial momentum p = (3 ■ E) 

{ 1 



'-total 



1 1 1 
- + - + - 

t U S 



iM 



tree 



1 1 1 
- + - + - 

t u s 



(23) 



This formula eventually coincides with what one gets from Eqns. (^JT^). 
The tree level amplitude is 



(24) 



and explicit use of Eqn. ([T^) allows us to arrive at the same Eqn. (pSf). 

As an alternative way of calculation, we computed /(4a) by a dispersive 
calculation. In this calculation we used the usual Landau rules to show 
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(a) (b) 



Figure 6: More infrared divergent one-loop graphs. 

that this graph has no other singularities at fixed t < in the complex s 
plane but a two- if cut from 4M^ to +00, corresponding to the cut graph 
on Fig. 0. The singularity for 4M^ < s < 4M^ — t is there in spite of the 
fact that the physical region for this process (we have fixed t < 0\) starts at 
s > 4M^ — t> 4M^ . The discontinuity across this cut is calculated from the 
Cutkosky rules 

/m{.A<,M} k . (25) 

A straightforward calculation leads to 

Im [iM^.a)} = -vr^ • 9' ■ j^TTT^ 7^ ■ log - + 0(1). (26) 

The real part of the amplitude is determined from the dispersion relation 

1 f°° ds' ( ~i 

iMibox) = -- -■Im\iM(^ia){s\t)\. (27) 

TT J4M2 S' — S — le ^ 

We note that no subtraction is needed and that we must include the 
contribution from unphysical s' < 4M^ — t where p is imaginary. Substitut- 
ing Eqn. (|26|) into Eqn. (|27|) some elementary complex integrations lead to 
Eqn. ( pTTl) . 



13 



cut 
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m 



m 



-• • — 

4M ' 4M '-t 



M 



M 



M 



(a) 



(b) 



Figure 7: Singularities of the graph (^a) in the complex s plane for fixed 



This calculation, especially the coincidence of Eqns. (|6| JT^) and (^31) shows, 
through a particular example, the correctness of what we did in Sect. |l[ 

IV Conclusion 

We have shown, both by general argument and explicit example, that in di- 
agrams involving scalar loops no powerlike infrared singularities occur. The 
introduction of spin does not alter this conclusion, and we use this result to 
show that the power coupling proof of the Equivalence Theorem put forward 
in 0] does not suffer from IR divergences. However, our reasoning allows 
us to arrive at this conclusion only in the massive sector, because the above 
analysis can break down when massless propagators are included into loop in- 
tegrations (for our purposes, gauge bosons and/or light Higgses of the MSSM 
model are considered 'massless'.) 

On the basis of these results we put forward the conjecture that, at least 
when the light particles are not coupled to each other — no worse IR di- 
vergences occur than O (^g"- ■ log" we actually prove this in Sect. y. A 
general proof seems straightforward. The generalization for coupled light par- 



t < 0. 
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tides (which is certainly the case in the SM) needs a separate investigation. 
This case is in a sense similar to the analysis of IR divergences in QCD. One 
particular graph that could cause trouble is the same as the one one Fig. |I], 
with light particles on the internal lines. Yang-Mills theories certainly have 
bad IR behavior which ultimately leads to confinement; the corresponing 
investigation might involve an adaptation the methods used there. 



Appendix: The calculation of 7(4^) 

In this Appendix, as an illustration to the main ideas in all these similar 
calculations, we briefiy describe the calculation of the IR divergent part in 
the box graph of Fig. |^a. The usual Feynman parameter integral is 

cx+y+z<i dxdydz 





2 4 

-IT ■ g ■ 



with 



w 



{Jm{x,y,z) -ie}"^ 

Jjn{x,y,z) = ■ {x + y) + ■ {1 - X - y)^ 
— s ■ z{l ~ X — y) — t ■ xy . 

A change of integration variables to u = 1 — x — y, v = 
= 1 — — allows us to have a form 

l-x~y 



(28) 



(29) 



^ and 



m 



[l-u) + M^-e + 



u 



(x+y) 



W 



V — S ■ U 



(30) 



in which all the u and v singularities are on the upper complex half plane. 
We change both contours to ones in the lower half plane as shown on Fig. ^ 
We then treat each of the resulting six integrals in turn. As an example, 
we take the parts with v = —iv : {—i 0) and w = 1 — iw : (1 — i — > 1): 

71-2 H rl 

iM. =^ ■ g"^ ■ du ■ u ■ (1 — u) ■ 

2 Jo Jo 



dv 



with 



■ {l-u)+u^ ■ (M^ 

2 



tv 



dw 



{Jm - ief 



(31) 



(32) 



— I 



u 
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(in v) 



(in w) 



Figure 8: Avoiding singularities in the Feynman parameters v and w on the 
complex plane. The dashed line represents the possible positions of singular- 
ities. 



(wc set the unit of mass to y/—t ^ !)• 

Next we break up the u integration as Jq^'^ du + du and (using Le- 
besgue's theorem) show that the second is IR regular. In the first term, a 
second order Taylor formula for allows to drop the remainder term and 

we are left with (dropping all IR finite terms) 

lAi =^ ' 9^ ' du ■ u ■ {1 — u) ■ I dw 

2 ^0 

.1 • (l + p) 

/o 



(33) 



Jo 



Jm{v = 0) 



1 2 



IR finite terms. 



An elementary decomposition of the integrand helps us to get rid of the 
square in the denominator 

f-i/^ u ■ du 



iM 



TT 



u 



dw - I 

Jo Jo 



dv 



J^{v = 0)-ie-i-(^ 
16 



(34) 
(35) 



+IR finite terms. 
We can now do two of tlie integrals to get 



iM 



ITT 



9 



1/2 

u 



du 



u 



2^2 " Jo 

[log (1 + z - ^) - log (1 + z + ^)] - [log I 



(36) 

- log {i + ^)] 



nil . _ _ ,V 

£2 „2 P 'fc 



Similar calculations for the region v : (— z — >■ 0) and w : (0 
to a similar formula with the numerator replaced by 



1 — lead 



log(v^+l + «) - log 



and all other regions end up with no IR divergent contributions. The sum of 
all IR divergent terms, with the variable x = is 



2^2 



dx f 1 
2 2x \ VI + 4x^ 



(37) 



log {i + log {i- + i 

^ [i^x — (3'^ — ie 



An analysis of the complex x singularities and changing the x contour 
shows that the ^^j^^^ term is IR regular. In the other term we use y = ^ ■ x 
and with similar tricks we can show that in the resulting formula, 



iM 



in 



AEp 



dy 



2^ y 



(38) 



log (^ + - 1 - i^) - log (^ - - 1 - i^) + 



^/y^^^^^ 

all the IR divergences are contained in a small neighborhood of zero, in an 
JJ' dy region with any m-independent r] > 0. Because the fraction in the 
integrand is an analytic function in < y < r], the IR divergence is read off 
easily: 



iM 



71 



AEp 



■9 



log 



l + (3 
1-/5 



/i2 



(39) 



which, restoring the unit of mass^, — 1 ^ t, yields Eqn. (|21.1| ). 
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